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ABSTRACT 

We model the mineralogy and distribution of dust around the white dwarf 
G29-39 using the infrared spectrum from 1-35 /im. The spectral model for G29- 
38 dust combines a wide range of materials based on spectral studies of comets 
and debris disks. In order of their contribution to the mid-infrared emission, 
the most abundant minerals around G29-38 are amorphous carbon (A < 8 /mi), 
amorphous and crystalline silicates (5-40 /mi), water ice (10-15 and 23-35 /im), 
and metal sulfides (18-28 /im). The amorphous C can be equivalently replaced 
by other materials (like metallic Fe) with featureless infrared spectra. The best- 
fitting crystalline silicate is Fe-rich pyroxene. In order to absorb enough starlight 
to power the observed emission, the disk must either be much thinner than the 
stellar radius (so that it can be heated from above and below) or it must have an 
opening angle wider than 2°. A 'moderately optically thick' torus model fits well 
if the dust extends inward to 50 times the white dwarf radius, all grains hotter 
than 1100 K are vaporized, the optical depth from the star through the disk is 
T|| = 5, and the radial density profile cx r -2 ' 7 ; the total mass of this model disk 
is 2 x 10 19 g. A physically thin (less than the white dwarf radius) and optically 
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thick disk can contribute to the near-infrared continuum only; such a disk can- 
not explain the longer-wavelength continuum or strong emission features. The 
combination of a physically thin, optically-thick inner disk and an outer, phys- 
ically thick and moderately optically thin cloud or disk produces a reasonably 
good fit to the spectrum and requires only silicates in the outer cloud. We dis- 
cuss the mineralogical results in comparison to planetary materials. The silicate 
composition contains minerals found from cometary spectra and meteorites, but 
Fe-rich pyroxene is more abundant than enstatite (Mg-rich pyroxene) or forsterite 
(Mg-rich olivine) in G29-38 dust, in contrast to what is found in most comet or 
meteorite mineralogies. Enstatite meteorites may be the most similar solar sys- 
tem materials to the G29-38 dust. Finally, we suggest the surviving core of a 
'hot jupiter' as an alternative (neither cometary nor asteroidal) origin for the 
debris, though further theoretical work is needed to determine if this hypothesis 
is viable. 

Subject headings: white dwarfs; stars: individual (G29-38, WD 2326+049); in- 
frared: stars 



Introduction 



The end state of stellar evolution for most ty pes of star is a white dwarf. Plane- 



tary material is common a round main sequence stars (IBackman fc Parescdll993l : iRhee et al. 



20071 : iTrilling et al.ll2008l ). Part of the planetary system is e xpected to survive the vio- 
lent, late stages of stellar evolution (IDebes fc Sigurdssonll2002l ). Thus it is to be expected 
that planetary materials are present around white dwarfs. But since white dwarfs are so 
under-luminous, this material is normally impossible to detect via reflected light or thermal 
emission, if it is distributed like the Solar System (with its inner portion destroyed out to 
at least the maximum radius of the red giant /AGB photospher e). Young, hot white dwarfs, 
where debris can be detected even at 100 AU from the star (ISu et al.l 120071 ). are the ex- 
ception to this rule. Planetary material can become detectable if it moves sufficiently close 
to a star, for example via gravitational perturbations such as produce comets in the inner 
Solar System (from their Kuiper Belt or Port c l oud h omes) and meteorites on the surface 
of Earth (from their asteroid belt homes). lJural (120031 ) explained how an asteroid can pass 
sufficiently close to a white dwarf to be tidally disrupted, leaving a disk of dust near the star. 
Sun-grazing comets in the solar system pass within similar distances. The new and growing 
class of dusty, metal-rich white dwarfs (classified as DAZd, where D=white dwarf, A=mostly 
hydrogen atmosphere, Z=trace atmospheric metals, d=dusty) allow a unique window into 
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planetary systems (Ivon Hippel et al.l 120071 ). 



The first-known DAZd star, and the one with the brightest infrared excess, is Giclas 
29-38 (WD 2326+049; G29-38 hereafter). Its effective temperature has been estimated 
in the range 11,800 to 12,100 K with su rfac e gravity logq in the range 8 .14 to 7.90 (e.g. 
Bergeron et al.ll2004l ; iKoester et al.ll2005l ). In iHolberg and Bergeron! (120061 ). the model cor- 
responding to the lowest estimated T highest g has radius 8100 km and mass .69 M , while 
the model corresponding to the highest T and lowest g has radius 9650 km and mass .55 
Mq. At wavelengths longer than 2.5 /im, the spectrum is completely domin ated by non 



photos pheric thermal emission with color temperatures 300-900 K. In paper 1 (jReach et al. 



2005al ). we measured the spectrum of this star from 5.5-14 /im, revealing a remarkably strong 
emission feature from 9-11 /im that is typical of small (< 2/ahi radius) silicate grains. The 
brightness of the disk permits spectroscopy at sufficient signal-to-noise to study the dust 
mineralogy. We have now measured the infrared spectrum out to 35 //m, and in this paper 
we discuss dust cloud models and compositions in more detail. 

After presenting the new observations, analyzing the variability, and compiling the spec- 
tral energy distribution in §2J we model the disk several different ways. In §31 we use a wide 
mix of minerals, and a simple temperature distribution, to determine likely constituents of 
the disk based on spectral features. In §H we present optically thin shell, moderately op- 
tically thick disk, and optically thick thin-disk models. In §51 we present constraints based 
on variability. In §H1 we tie together the evidence gleaned from comparing the models to the 
observations. 



Observations 



Table [T] shows a log of observation s with the Spitzer Space Telescope ( IWerner et al. 



20041 ). The Infrared Spectrograph (IRS) (IHouck et al.l 12004 ) observations were performed in 
"staring" mode, wherein the source was placed on each of two nod positions on each of the 
slits and orders. IRS data were processed with the S16 pipeline (with backgrounds subtracted 
by differencing nods) and extracted using optimal extraction in the Spitzer Science Center's 
data analysis tool SPICE. The 4.5, 8, and 24 /im observat ions described in Paper 1 were 
supplemented by archival Infrared Array Camera (IRAC) (IFazio et al.ll2004l ) 3.6, 4.5, 5.8, 
and 8 /im observations. Annular aperture photometry was performed on each IRAC basic 
calibrated image with a 4 pixel on-source radius and 8-20 pixel background annulus; array- 
location-dependent, aperture- loss, and pixel-phase photometric corrections were applied; and 
the uncertainty of each image's photometry was d etermined by comb ining uncertainties due 
to photon statistics and background removal (see lReach et al.ll2005bl ). 
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2.1. Variability of the mid-infrared flux 

The variability of mid-infrared flux of G29-38 can be constrained with the IRAC photom- 
etry. Fluxes are shown in Table [2j The brightness of the star on each basic calibrated image 
from each observing sequence was measured, corrected for the array-location-dependent re- 
sponse, then a weighted mean and statistical uncertainty in weighted mean computed. The 
absolute calibration uncertainties are not included, because we are comparing fluxes from 
the same instrumen t; IRAC photometry has been shown to be stable to better than 1% 



( iReach et al.l l2005bl ). To obtain the highest signal-to- noise from each image for repeata- 
bility, we found an aperture radius of 4 pixels was optimal. The signal-to-noise was 230, 
2008, 90, and 130 at 3.6, 4.5, 5.8, and 8 /mi. At 4.5 & 8 jttm, the star was observed in two 
independent observations just a few minutes apart, yielding identical mean fluxes; the flux 
differences are —0.1 ± 0.2% and —0.1 ± 0.3%, respectively. At 3.6 & 5.8 /mi, the star was 
observed with different AORs separated by 13 months. The star was fainter in the second 
epoch, by —3.3 ± 0.4% at 3.6 /im (and a statistically insignificant —1.0 ± 0.8% at 5.8 /im). 
The brightness difference at 3.6 /mi is likely due to short-term fluctuations in the star, to 
which we turn now. 

The variability of the mid-infrared flux over shorter time scales can be asses sed from the 



individual images taken in 2004 Nov. By the design of IRAC (IFazio et al.ll2004j ). the 3.6 and 
5.8 /tin channels observe simultaneously, as do the 4.5 and 8 /tm channels. The sequence of 
events was as follows, with times given relative to the first image at 2004-Nov-26 10:54:11.8 
UT: 



00:00 the star is placed in the 4.5+8 /im field of view and a short (1.2 sec) frame is taken, 

00:03 5 consecutive dithers, consisting of a 30 sec frame and a short telescope slew, are 
performed, 

04:34 the telescope is moved to re-center the star in the 4.5+8 /tm field of view and a short 
(1.2 sec) frame is taken, 

04:37 19 consecutive dithers, consisting of a 30 sec frame and a short telescope slew, are 
performed, 

18:22 the telescope is moved to center the star in the 3.6+5.8 /tm field of view, 



1 We tested the IRS observations for variability, but there was no clear evidence. The IRS observing 
strategy provided 6 samples along 60 sec ramps in each subslit. The flux varied by < 4% among the shortest 
wavelength spectra, with no wavelength dependence detectable above the noise. 
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18:22 19 consecutive dithers, consisting of a 30 sec frame and a short telescope slew, are 
performed, then 

30:51 the sequence is complete. 

Figure [1] shows the flux versus time. 

The scatter in the photometric measurements is significantly larger than their uncertain- 
ties, especially at 3.6 /on. To ass ess the harmon ic content of the time series, we computed 



the periodogram as described by IScargld (119821 ). Figure [2] shows the periodograms at the 



4 IRAC wavelengths. The time series cannot be described by a simple period, but instead 
contain variation on a range of timescales. This is characteristic of ZZ Ceti variable stars 



(IKleinman et al.lll998l ). The peak in the 3.6 /im amplitude at 190 sec is highly significant 
(14a), as is the harmonic content at 300-440 sec. None of the other harmonic content is 
significant, at the sampling frequency and signal-to-noise level of these observations, except 
some power (3.3a) at ~ 220 sec at 4.5 /im. 



The flux variations seen with IRAC are similar to those seen at Palomar by lGraham et al. 



(119901), who observed simultaneously at wavelengths dominated by the photosphere (B, J) 
and the infrared excess (K-band). They found fluctuations at periods 181 and 243 sec in 
K-band, with no corresponding ones at J-band. Our IRAC 3.6 /im results confirm the sig- 
nificant fluctuations at ~ 190 sec. The amplitude of the fluctuations ~ 4% at 3.6 /xm, while 



Graham et al.l (119901 ) found 2.5% variations at 2.2 /xm. The stellar photosphere must be 
subtracted before interpreting these results. At 2.2 (3.6) /im, the photosphere contributes 
66% (18%) of the total flux. If we assume the infrared photosphere is constant, and the 
fluctuations are due to the disk, then the amplitude of disk fluctuations at 2.2 (3.6) fim is 
7% (5%), i.e. very similar, with K-band possibly slightly higher in amplitude. Fluctuations 
at 4.5, 5.8, 8 /im are not detected in the periodogram, with upper limit ~ 5%, 5%, 3% of 
the total flux, or 5%, 5%, 3% of the infrared excess above photosphere (2a limits). The 
upper limit at 8 /im, and the trend such that fluctuations are most significant at 2.2 and 3.6 
/im, with a rapid decrease in amplitude at longer wavelengths, are significant and indicate 
that the fluctuations are due to an emitting region with a relatively high (> 1000 K) color 
temperature. 



2.2. Compilation of the spectral energy distribution 

In order to model the mineralogy and dust distribution around G29-38, we need to 
combine the observed data from the near- through mid-infrared, correcting for calibration 
errors as well as source variability. This is in fact not possible, since the temporal sampling of 
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the photometry is inadequate, and the observations are not contemporaneous. We therefore 
must assemble the various portions of the spectrum in such a way as to make them most 
plausibly 'connected'. In principle, there is a scale factor for each wavelength range that 
depends on the epoch. The amplitude of the scale factor, due to source variability, may 
reach up to ~ 10% in the near-infrared and should be less than 5% at longer wavelengths. 
Calibration uncertainties range from ~ 3% in the optical to 5% for IRAC and MIPS to 10% 
for IRS. The range of plausible scale factors is then 5-10% across all wavelengths. 

We first ensured that the spectroscopy and photometry were in accord. The IRAC fluxes 
must be corrected for the spectral shape of the source first. For the IRAC 8 fim channel, 
which is fully cover ed by the IRS spectrum, we integrated the IRS spectrum appropriately 



( IReach et al.ll2005bl ) over the passband to calculate the color-correction K — 1.16 (where the 
corrected flux is the 'observed' flux from the IRAC calibration divided by K) and applied this 
correction to the photometry. The IRAC 8/im flux density can then be compared directly 
to the IRS spectrum: the flux at the nominal wavelength 7.872 /im is 7.19 mJy from IRAC 
and 7.30 mJy from IRS, a deviation of only 1.5%, well within the calibration error budget. 
The IRAC 3.6, 4.5, and 5.8 jum flux density required only very small color corrections: the 
spectrum has a color temperature ~ 1100 K at these wavelengths, for which the corrections 
are 0.995, 0.996, and 0.996, respectively. The absolute fluxes from IRAC and IRS at 5.831 /im 
are 7.76 and 7.54 mJy, showing the instrument cross-calibration is good. The observed MIPS 
24 fj,m flux was divided by a color correction of 0.97 to account for the wide spectral response 
of the filter and detector. For IRAC and IRS, we therefore find that there is no need for a 
relative re-scaling; that is, the measured flux densities averaged over their exposure times at 
their epoch of observation agree well, whether due to averaging of or lack of variability. 

To remove the white dwarf photosphere from the spectra, we use a model atmosphere 
for T = 12, 000 K and log g = 8 covering 0.35-60 /im (courtesy D. Koester), normalized using 
2MASS photometry (J = 13.132 ± 0.029, H = 13.075 ± 0.022 K, = 12.689 ± 0.029) and 



optical spectrophotometry from Palomar (IGreenstein fc Liebertl Il990l ) . In fact, the optical 



and 2MASS photometry cannot both agree with the model spectrum. If we normalize the 
white dwarf model at the 2MASS J band, then the model would unde r-predict the visible 



flux by 8%. The near-infrared spectrum obtained by IKilic et al.l ( 120061 ) with the IRTF also 
disagrees with the 2MASS photometry; if we normalize the white dwarf model by the 2MASS 
J-band flux, then we under-predict the IRTF spectrum by 18%. The low flux observed by 
2MASS is likely due to variability of the white dwarf. The shape of the IRTF spectrum 
follows the model closely at wavelengths shorter than 1.6 /im; the infrared excess is evident 
at wavelengths beyond 1.7 /im. To obtain a joined spectral energy distribution we proceed as 
follows: (1) normalize the white dwarf model to match the optical photometry, (2) scale the 
2MASS photometry by a factor 1.087 to match the model at J band, (3) rescale the IRTF 
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spectrum by 0.92 to match the white dwarf model at 0.8-1.3 /xm, (4) rescale the 2MASS 
photometry and IRTF spectra by a factor of 1.2 so that the red end of the IRTF spectrum 
(2.4 /mi) is plausibly consistent with the IRAC 3.6 /xm photometry. (Plausible consistency 
here was defined as allowing a blackbody to pass through the IRAC 3.6 /im and 2MASS 
K s photometry for color temperatures 800-1200 to within the error bars, and keeping the 
IRTF and 2MASS photometry consistent at 2.17 //m.) The uncertainties were taken to be 
a root-sum-square combination of the measurement uncertainties and 10% of the brightness 
of the photosphere, to account for scaling uncertainties in the photospheric subtraction. 

Figure [3] shows the combined spectrum of the infrared emission after subtracting the 
stellar photospheric emission. The combined scalings yield near-infrared photometry and 
spectrum effectively at the same epoch (in terms of variability) as the mid-infrared photom- 
etry. 



3. Mineralogy 



The composition of the dust around G29-38 was determined by comparing the spectra 
to those of a set of constituent materials, using a fitting method that has been used for comet 
Hale-Bopp, the Deep Impact ejecta from c omet 9 P /Temp el 1, and th e dust around the stars 
HD 100546, HD 69 380, and HP 113 766 Lisse et all liooS . l2007al lbl. boosh . Details of this 
model are given in iLisse et al.l (120081 ). In this section, we model the emitting region as an 
optically thin dust torus. We apply this simple model to G29-38 to allow direct comparison 
to dust in other astronomical systems, but we will revisit the mineralogy below when testing 
more sophisticated models. 



3.1. Particles sizes and temperatures 



A toroidal mod el is motivated by the s ingle temperature distr ibutions found for many of 
the stars studied by lBeichman et al.l (120061 ) and I Chen et al.l (12006!) as well as the narrow dust 
structures found in many of the HST images of debris disks (IKalas et al.ll2005l 120061 ). The 
best-fit temperature of the smallest dust particles (0.1-1 /im), which superheat substantially 
above Local Thermal Equilibrium (LTE), is T max = 890 K for the olivines, 850 - 890 K for the 
pyroxenes, and 930 K for the amorphous carbon. The largest particles in our calculation, 
with radius 1000 /im, are set to T^tr = 6 00 K, and dust of intermediate size is scaled 
between the two extremes (ILisse et al.l 120061 ) . Using the Tempel 1 ejecta temperatures as 
a guide (where we found olivines at 340 K and amorphous carbon at 390 K at 1.51 AU 
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from the Sun, where Tlte = 230 K), and allowing for an G29-38 stellar luminosity that 
is 2 x 10~ 3 L Q , we estimate a location for the hot dust ~ 150 stellar radii from the white 
dwarf. The location of colder dust, capable of supporting a stable water ice component at 
200 K, is at > 10 3 stellar radii (> 8 x 10 11 cm). While differing by at least a factor of 9, 
this range of distances is still small, suggesting material in tight orbit around the WD. The 
best-fit single continuum temperature to the 7-35 y^m Spitzer spectrum is 950 K, close to 
the amorpho us carbon tempe rature. As seen in the Tempel 1 ejecta from the Deep Impact 



experiment (ILisse et al.ll2006l ). the amorphous carbon dominates the continuum behavior 
as it has a featureless emissivity and is the hottest material, thus contributing most to the 
short-wavelength emission. 

The best-fit size distribution dn/da oc a - 3 - 7±0 - 2 argues for dust surface area dominated 
by small particles, but dust mass dominated by large particles. A system in collisional equi- 



ff 

librium would demonstrate a PSD oc a 3 5 (jDohnanyj Il969l ; iDurda and Dermottl 119971 1; for 



"real" systems a size distribution even steeper than a at small si zes is expected in a colli 



sional cascade because of the dependence of particle strength on size (jO'Brien and Greenberg 



2003|). 



The total dust mass required to explain the Spitzer IRS spectrum (i.e., mass in particles 
of 0.1-10 /im in size that contribute appreciably to the xl value of the fit to the infrared 
spectrum) is ~ 2 x 10 19 g. Extrapolating (using the best-fit size distribution) from a maxi- 
mum size of 10 /im to 1 cm would increase the mass by a factor of 8. A total cloud mass of 
order 10 19 g compressed into a solid sphere of average density 2.5 g cm -3 (i.e. rocky silicate 
material) would have a radius of 10 km, equivalent to a single, small asteroid or large comet. 
The surface area of detected particles is 5 x 10 22 cm 2 . If this dust is in an annulus of inner 
radius 7 x 10 10 cm (90i?*) extending to twice that radius, then the areal filling factor of 
grains viewed from above the disk is of order unity. The optically-thin model is thus unlikely 
to apply, although most of the observed emission (in particular the spectral features) must 
arise from the optically thin parts of the disk. We address the optical depth effects in the 
modeling sections below, but we proceed first (with caution) to discuss the mineralogy from 
the optically thin fits. 



3.2. Dust Composition 



Over 80 different species were tested for their presence in the spectra. The material 
spectra were selected b y their reported p resence in meteorites, in situ comet measurements, 
YSOs, and debris disks ( Lisse et aDbood ) . Consultations with members of the Stard ust team, 
and examination of the interplanetary dust particle (IDP) literature (cf. review by iBradley 
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20021 ). and the astrominero logical literature (cf. review by lMolster and Wateral2003l ) pointed 
to the most likely mineralogical candidates to be found in the dust. Materials with emissions 
matching the strong features in the Spitzer emissivity spectra were also tested. The list of 
materials tested against the Spitzer spectra included multiple silicates in the olivine and py- 
roxene class (forsterite, fayalite, clino- and orth-enstatite, augite, anorthite, bronzite, diop- 
side, and ferrosilite); phyllosilicates (saponite, serpentine, smectite, montmorillonite, and 
chlorite); sulfates (gypsum, ferrosulfate, and magnesium sulfate); oxides (various aluminas, 
spinels, hibonite, magnetite, and hematite); Mg/Fe sulfides (pyrrohtite, troilite, pyrite, and 
niningerite); carbonate minerals (calcite, aragonite, dolomite, magnesite, and siderite); wa- 
ter ice, clean and with carbon dioxide, carbon monoxide, methane, and ammonia clathrates; 
carbon dioxide ice; graphitic and amorphous carbon; and polycyclic aromatic hydrocarbons 
(PAH). Of these materials, a small, unique subset was found necessary to properly fit the 
Spitzer data. Sources for the data included: for s ilicates, the Jena Datab ase of Optical Con- 
stants for Cosmic DustH ( jDorschner et al. 1995 ; Jager et al. 1998 . 2003 ). the Mars Global 
Surveyor Thermal Emission Spectromete r databasqj, and W . Gl accum (private comm uni- 



cation), as well as emission spec tra from iKoi 



carbonates, Kemper et al.l (120021); f or sulfides 



Keller et al. 



for amorphous carbon, lEdohl (119831 ); for PAH, iLi and Draind ( 120011 ) 



ke et all (120031) andlChihara et al.l (120021); fo r 



(2002) and 



Kimura et al 



(120051 ): 



We determined the reduced x 2 f° r the model fit to the data for thousands of combinations 
of minerals. All models with reduced chi-squared values large than the 95% confidence limit 
were excluded from consideration. The range of constituents abundances was determined by 
varying the amount of a material, and finding where the model exceeded the 95% confidence 
limit. Only constituents with abundances significantly above zero were included in the final 
best fit. Table [3] shows the composition of the best-fit model. For each entry, we also show 
the reduced \ 2 f° r the model if that constituent is deleted. Upper limits for non-detected 
species are included in Table Abased on the 95% confidence level amplitude of the constituent 
when it is included in the fit. 

Figure 0] shows the best-fit spectral model for G29-38. In order of their contribution 
to the mid-infrared emission, the contributing minerals are amorphous carbon; amorphous 
olivine; crystalline silicates ferrosilite (FeSi0 3 ), fayalite (Fe 2 Si0 4 ), diopside (MgCaSi 2 6 ), 
and enstatite (MgSi0 3 ); and metal sulfides (Mg 10 Fe 90 S); and water ice. Figure \5\ shows the 
spectral model after removal of the best-fitting amorphous silicates that dominate the overall 
spectral shape. The fit to the IRS data (5.2-35 jttm) is excellent, with xl — 1-03. 



2 http: / /www. astro. uni-jena.de/Laboratory/OCDB 
3 http://tes. asu.edu 
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The composition of the dust surrounding G29-38, as determined by modeling of the 
Spitzer IRS spectrum, is unusual when compared to circumstellar dust in other environments. 
Amorphous olivine is present, as in the other systems that have been modeled with the same 
technique. But the relative lack of crystalline olivine, and the strong prevalence of Fe-rich 
silicates, make the G29-38 material distinct from other systems. The mix of crystalline 
pyroxenes and amorphous olivines may be indicative of 'aged' dusty material. No PAHs, 
carbonates, or water gas are seen, reflecting a total lack of primitive nebular material. Some 
metal sulfides appear to be present, suggesting temperatures as low as ~600 K in the observed 
dust; 600-700 K is the temperature range for vaporization/condensation of ferrosulfides. 
Given that the best-fit continuum temperature for the spectrum is 930K (dominated by the 
short wavelength emission from 0.1 - 1.0 /xm amorphous carbon particles), either the carbon 
is superheated beyond the temperature of the metal sulfides, which may be present in larger, 
cooler particles, or there is a distribution of dust locations and effective temperatures, i.e. a 
disk-like structure. 



3.3. Water ice 

The detection of water ice emission in the spectrum is curious. In a vacuum, water 
ice sublimates at temperatures above ~200 K. Thus water ice cannot be in direct physical 
contact with the rest of the hot dust detected in the infrared spectrum. Either the dust 
must be continuously created, or it must reside in a physical location removed from the rest 
of the dust reservoir. Since no water vapor is detected at 6 fim, where there are strong 



features that should be dete cted if water were present in significant quantities (ILisse et al. 



20061 ; IWoodward et al.l 120071 ). no appreciable ongoing sublimation can be occurring. Water 
vapor is i onized within ~ 10 3 sec, if w e scale from the lifetime of water at 1 AU from the Sun 
(~ 10 6 s; ISchleicher fc A'HearrJfl98~8l ) according to the luminosity of G29-38 (2 x 1O~ 3 L ). 



There would still be a steady state amount of water vapor present if the ice is sublimating. 
Thus it appears that the water ice is ctt £i lcir ge remove from the rest of the warm circumstellar 
dust, and the water ice is found in the IRS beam (i.e. within 5 x 10 14 cm [7 x 10 5 .R*] of the 
star). 

In paper 1, we showed that a range of temperatures is required to match the photometry 
out to 24 /iin, with a two-temperature fit having 890 and 290 K color temperatures. Thus we 
already suspected either a continuous range of dust temperatures or the presence of two dust 
reservoirs, one hot and one cold. For G29-38, the ice (or ice-coated grains) must be located 
at > 10 12 cm (2000-R*)from the star, which, while farther than the dust that dominates the 
mid-infrared spectrum, is still very close to the star. Cold, icy dust around other stars, at 
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locations comparabl e to the Solar Syste m's Kuiper Belt, is common, being present around 
15-20% of all stars (IBryden et al.ll2006l ). The ice detected around G29-38 is closer to the 
star than the far- infr ared debris disks at ~ 100 AU (10 15 cm) radius commonly seen around 
main sequence stars (IBryden et al.ll2006l ). and the present location of the icy grains would 
have been within the star when G29-38 was in its red giant or AGB stage. The source of 
this cold, presumably icy, material is likely to be the same parent bodies that create the 
dust closer to the star. We will return to the possible parent body natures in the Discussion 
section, but we should point out now that presence of H 2 is well established in asteroids 
(especially outer belt), meteorites (as water of hydration as well as signatures of aqueously 
altered mineralogy), comets, and planets. The H 2 in any of these bodies would freeze if it 
were liberated by disintegration of these bodies and survive out to > 10 12 cm from the white 
dwarf, owing to its small luminosity. 



4. Distribution of material around the star 

Two possible disk models are considered; these are illustrated in Figure [6] and discussed 
in turn in the following sections. The models will be fitted to a combination of the Spitzer 
spectroscopy and photometry, the 2MASS photometry, and the near-infrared spectrum, after 
removal of the photosphere model and scaling to a common epoch (§2.2). 



4.1. Physically thick disk 

If the optical depth through the dust cloud around G29-38 is optically thin at the 
wavelengths of the observed emission, then the spectrum is determined by simply integrating 
the density distribution, weighted by the Planck function at the local temperature, through 
the cloud. If the cloud is optically thin at visible wavelengths, where the spectrum of the 
white dwarf peaks, then the dust heating is simply determined by the integrated absorption 
of the star's distance-diluted spectrum by each grain. We can place some constraints on 
whether such a model can apply. First, the cloud can only emit as much energy as it absorbs 
from the star. Let us consider a flattened torus, defined in spherical coordinates r, 9 (with r 
distance from the star, 9 the angular separation from the equatorial plane) as having non- 
zero density for R\ < r < R 2 and 9 < 9i. Since the dust luminosity is / = 3% of the 
star luminosity, the opening half-angle of the torus must be at least 9i > f/2 radians, i.e. 
9i > 0.8°. A thinner torus simply cannot intercept enough starlight to emit the observed 
luminosity. (The constraint does not apply to a disk thinner than the radius of the star, nor 
to a warped disk, as discussed in following sections.) 
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Of critical importance to the radiative transfer for calculating the emergent spectrum 
from the disk is its optical depth in the infrared. For a torus with radial mass density profile 
p = pi(r/i? 1 ) _Q? , the total cloud mass 

M = 4n9 hPl Rlf a , (1) 

where 

(R 2 /Ri) 3 ~ a - 1 , , , , 

/a = ^ (for«^3) (2) 

o — a 

= \n(R 2 /R 1 ) (fora = 3). (3) 

The vertical optical depth at the inner radius is 

3QirM 
±X l6irp d aRlf a 

where a is the particle radius and Qm is the the absorption efficiency averaged over the 
thermal emission spectrum. Using the cloud mass radius from paper 1 (Ri = 1R & ) and 
scaling the mass in units of 10 18 g, t_i_i = 0.03Mi$Qjf>(a/ /im)" 1 . For particles smaller than 
the wavelength, the absorption efficiency can be approximated as Qm — 27ra/A, so the 
optical depth for thermal emission at wavelength A is r^ 1A = 0.2M 18 (A//zm) _1 , independent 
of particle size. For thermal emission at 10 /im, one can neglect radiative transfer out of the 
disk only if M < 5 x 10 19 g. 

Dust heating is determined by the propagation of starlight through the disk. The optical 
depth from the disk inner boundary to a distance twice as far from the star is 

T ±iQopt 1 — 2 a 
11 ~ 6iQir a-1 • 1 ) 

For particles larger than 0.1 /mi, Q v — 1 for absorption of starlight, and the optical depth 
from the star to the disk interior is Tp = 2.5M 18 (9i /OS )^ 1 a~^(Ri / R & )~ 2 . Using the mass 
and radius from paper 1, we find the cloud is optically thin for starlight propagation (and for 
thermal emission) as long as Oia^ > 2°M 18 . The presence of the silicate emission feature 
further requires the particle size < 2/im, so for a disk that is optically thin to propagation 
to starlight, we require M 18 < 6i. The cloud must be at least as massive as derived in paper 
1 (where light was allowed to propagate through the disk unimpeded), so we finally obtain 
the constraint 1 < M 18 < 0i for an optically thin disk. If the disk is more massive than the 
upper limit, it becomes optically thick to propagation of starlight, the temperature decreases 
more rapidly with distance from the star, and the cloud is more dominated by its inner edge. 
We will refer to this case as moderately optically thick, with 0.5 < tii < 10 so the cloud is 
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optically thin to vertical propagation of the observed infrared emission but optically thick 
to radial propagation of starlight ( §4. 1 .2[) . 

It appears therefore that a model cloud that is 'physically thick' (i.e. thicker than 
the white dwarf's diameter) can still be quite flattened, with angular widths as small as 
1°, before becoming optically thick. The distinction between t he 'physically thick' model 
( jReach et al.ll2005al ) and the 'physically thin' model (1 Jural 120071 ) is that the physically thin 
model is thinner than the stellar diameter, so that it can be illuminated on its upper and 
lower surfaces. In order to intercept as much energy as is observed (3% of the white dwarf's 
luminosity), a physically thick model at its inner edge must have a half-width at least 10 9 
cm, which is larger than the stellar radius (~ 8 x 10 s cm). Starlight does not illuminate the 
surface of the physically thick model, and instead the dust is heated by stellar photons that 
must propagate through the disk. 



4-1.1. Spherical shell model 

The simplest geometry for the cloud around G29-38, and a limiting case for the physi- 
cally thick model, is a spherical shell. A moderately optically thick spherical shell will have 
the same temperature distribution as a disk with the same radial density variation. We 
use the spherical shell calculations to measure the radial temperature profile relative to the 
optically thin case. 

Spherical shell calculations were performed using DUSTY for a radial profile r 3 , an 
amorphous carbon or silicate composition, and a range of total optical depths, r (at 0.55 
/im), from 0.01 to 10. The inner boundary of the shell is where the dust temperature reaches 
1100 K. (Fits with inner temperature 1000 and 1300 yielded significantly worse fits to the 
observed spectrum.) Figure [7| shows emergent spectra compared to the observations. The 
first are visually very good, clearly reproducing the spectral shape, despite the simplicity of 
the model; however, the reduced xl — 2-2 is statistically poor. The residuals are dominated 
by structure near the 9-11 /im silicate feature: the observed minus model residuals have 
positive peaks at 9.2 /im and 11.2 /im and a broad negative trough spanning 8-13 /im. The 
deficiencies of this model are thus due to using only one silicate, in contrast to the more 
detailed mineralogy found in the multi-composition models discussed above. 



4 Ivezic, Z., Nenkova, M. & Elitzur, M., 1999, User Manual for DUSTY, University of Kentucky Internal 
Report, accessible at http://www.pa.uky.edu/~moshe/dusty 
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The luminosity of the shell per unit stellar luminosity is accurately approximated by 

^ = / ~ / (l - e -) (6) 

where /o=0.74 for carbon and 0.57 for silicate grains. The observed luminosity ratio, / = 
0.03, could be explained by the combined carbon and silicate cloud from Figure[7]with optical 
depth at 5500 A and 2000 A of 0.028 and 0.044, respectively. The total extinction toward the 
star is small but could potentially be measured with precise UV/visible spectrophotometry, 
in which case presence of carbonaceous grains could be test ed by searching for a 2175 A 



feature such as seen in the interstellar medium (IMathia Il994f ) . 



The temperature at a given distance from the star is the same for a spherical shell as 
for a moderately optically thick disk, as long as scattered light and dust thermal emission 
is a negligible heat source. Figure [S] shows the temperature profiles through shells with 
different optical depths. The temperature versus distance from the star for grey grains in 
the optically thin limit would follow T oc r~ 05 . For real materials, and taking into account 
radiative transfer, the temperature profiles are significantly different. Using a power- law 
approximation T oc r~ s , the predicted temperature profiles have 0.44 < S < 0.48 as long as 
r < 0.2. But for higher r a single power law is not adequate. We fitted the curves in Figure[H] 
with empirical functions, for use in the moderately optically thick models. For reference, 
the temperature scalings are closely related to t he parameter \l/ defined in the original paper 



on the self-similarity solution used in DUSTY (jlvezic and Elitzurl 119971 ). 



4-1.2. Moderately optically thick model 



Armed with the temperature profiles, and the wavelength-dependent cross-sections for 
each mineral, we can now compute the brightness of a flattened cloud that has optical depth 
7~|| < 10 and tjl < 1. We define this as a 'moderately optically thick' disk. For simplicity we 
consider a fanned disk, with scale height proportional to distance from the star, h = r tan 9 1 . 

2 

The model does not apply to a thin, flat disk or a disk with a vertical density gradient. Full 
treatment of the optically thick disk with vertical density gradient requires a two-dimensional 
calculation that is beyond our present scope. 

The moderately optica lly thick model wa s calculated for a sub set of the minerals use d 



in §3) amorphous carbon (jZubko et aL 



1996 ), amorph ous olivine (IDorschner et al.lll995l ). 



forsterite (pure-Mg crystal line olivine. I Jager et al.ll2003l ). and enstatit e (pure-Mg crystalline 
pyrox ene J Jager et al.lll998l ). and bronzite (Fe-rich crystalline pyroxene. iHenning and Mutschke 
19971 ). Based on the results of §31 we expect the main constituents (in order) to be the amor- 
phous silicates, carbon, and bronzite. Due to the moderately intensive computations and 
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lack of in-hand UV-FIR laboratory data (required to span absorption of the white dwarf 
spectrum as well as thermal emission), we did not include as wide a range of minerals as 
in §21 For bronzite, we use the calculations for forsterite but scaled them by the relative 
small-particle emissivity over the range of wavelengths where the optical constants were 
available (6.7-500 jttm). The most important compositions for which we did not calculate 
the moderately-optically-thick model were magnesium-iron sulfides, and water ice (which in 
fact cannot exist on grains in this model and requires a separate reservoir, as discussed in §3]). 
The temperature of grains of each size (from 0.1 to 1000 jttm) was calculated by balancing 
radiation from the white dwarf, geometrically diluted by r~ 2 within 20 < r/R* < 70000, with 
the grain's thermal emission. Figure [9] shows the grain temperatures for two materials and 
three particle sizes in this optically-thin limit. The temperatures were then adjusted using 
the scale factors derived in the previous section appropriate for the composition, distance 
from star, and total cloud optical depth. If a grain's temperature exceeds a vaporization 
temperature, T vap , it's emissivity is set to zero. Values of the vaporization temperature in 
the 1000-2000 K range are expected for most minerals; evidence for inner edges of YSO 
disks at distances corres ponding to these vaporization te mperatures has been found in inter- 



ferometric observations (IMonnier fc Millan-Gabetl 120021 ). The emission spectrum of dust at 



each distance from the star was calculated by integrating over several size distributions: a 
power-law oc a -3 5 , a Hanner law with slope 3.7 and critical size 7.4 /im, the size distribution 
from the coma of c omet Halley, and the size distribution of interplanetary meteoroids and 



lunar microcraters (jGriin et al.lll985l ). 



The emergent spectrum from the cloud was then calculated by integrating in spherical 
coordinates using a radial density distribution n oc r~ a and minimum radius Rmin-, with 
individual models sampled from the ranges 0.3 < a < 6 and 50 < R m i n /R* < 1000. We 
assume azimuthal symmetry and optically thin infrared emission, so the opening angle of the 
disk 8i does not affect the spectral shape. A range of 1000 < T vap < 1600 K and < t\\ < 10 
were considered. Then for each cloud geometry (a, R m i n ), a linear combination of the models 
for carbon, olivine, and forsterite was fitted to the observations, and the x 2 for the mixture 
computed. Because these materials vaporize at R < 507?* for essentially all particle sizes 
(see Fig. [9]), all models with R m in < 50-R* are equivalent. 

The best-fitting solution is enumerated in Table HI Given 390 data points in the spec- 
tral energy distribution, we expect a 'good' fit to have reduced xl — 1 within dispersion 
a/2/390 = 0.07. The best-fitting model had xl — 1-23. We discuss the residuals, which 
are localized rather than spanning a wide wavelength range, in the following paragraph; for 
statistical purposes we assume they are unrelated to the cloud geometry. The confidence 
intervals for the parameters, determined by Axl = 0.07, are listed in Table HI 
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The residuals from the radiative transfer model are shown in Figure [TTJ Coherent fea- 
tures dominate the residuals; the amount by which they increase xl measures their statistical 
significance. These include the following: (1) a 'W'-shaped pattern between 8 and 10.5 /im 
(Axl = 0.30); (2) an 'S'-shaped pattern between 15 and 21 /im (Axl = 0.24); and (3) a 
'U'-shaped residual between 1 and 3 /zm (Axl = 0.20). Each of these features is statistically 
significant. Feature (1) is due to mismatching the fundamental silicate Si-0 stretch feature; 
an additional silicate is required which provides a shorter-wavelength silicate feature than 
Mg-rich olivine. Feature (2) is due to either mismatching the silicate Si-O-Si bending feature 
or presence of MgO (or Mg x Fei_ x O with x > .5), which has a peak at 18.5 /im. Cosmic 
abundances favor formation of Mg-rich silicates over Mg oxides, so a preferable solution for 
residual features (1) and (2) is a silicate with bluer fundamental and better-matching Si-O-Si 
bending mode features. 

We experimented with the following minerals, inspired by the results of §31 All models 
included amorphous carbon and olivine. Additional (third) minerals were included one at 
a time; the most successful were bronzite (xl = 1.22) and forsterite (xl = 1.55), in order 
of goodness of fit. These results are consistent with those in §31 where an Fe-rich pyroxene 
was shown to be the third-most important mineral. Then for the most successful third 
minerals, fourth minerals were added; a good combination was forsterite plus montmorillonite 
(xl = 1.23), but it was only slightly better than bronzite model alone. Other minerals 
such as amorphous pyroxene yielded negligible improvement to the fit (Axl = 0.02). The 
materials commonly referred to as amorphous olivine and amorphous pyroxene are both 
amorphous silicates and are essentially indistinguishable. It appears to be significant that Fe- 
rich pyroxene is more abundant than Mg-rich olivine; this is directly due to Fe-rich pyroxene 
having a bluer fundamental Si-0 stretch (as observed) compare d to forsterite. From the 



laboratory measurements of phyllosilicates by lGlotch et al.l (120071 ). there are two properties 



that make them possible explanations for residual features (1) and (2): phyllosilicates have 
an Si-0 stretch that is bluer than forsterite, and they have a double-peaked Si-O-Si bending 
mode that is similar to the shape observed in the G29-38 spectrum from 18-20 /im. 

The distinction between the bronzite model and the forsterite+montmorillonite model 
is mathematically within the errors, but we will refer primarily to the former model since it 
is simpler and it agrees with the more extensive mineral search from §31 The fit could be 
improved if a more exhaustive mineral search were performed. However, the results of §3] 
already demonstrate that including a wide range of minerals can explain the silicate feature 
shape adequately, and we are already at a reduced xl — 1-2 so that the present signal-to-noise 
does not allow further, unique modeling. 

Feature (3) in the residuals relates to the shape of the inner edge of the disk and the 



-17- 



vaporization temperature. We only explored T vap on a grid with 100 K intervals, and we 
fixed it at the same value for all compositions, so the observed, modest deviations from our 
model are not particularly surprising in the near-infrared. Further, we did not explore the 
composition of the featureless materials, which could affect not only their T vap but also the 
slope of their absorption around 3-6 /mi, which determines the shape of the inner edge of 
the disk spectrum. 



4.2. Physically thin disk 



A physically thin disk is very optically thick, so starlight cannot propagate radially 
through the disk. Such a disk must be so thin that the star can illuminate its surface, or it 
must be warped such that the surfa c e has clear lines of sight to the star, or a combination 



Jura 



of both effects as described by 
scales as r~ 75 (IChiang fc Goldreich 



(2007). The temperature versus distance from the star 
19971 . e.g.). The flux from an optically thick disk is 
straightforward to estimate: for G29-38, the disk temperature T = 8008(-R/-R*)~ a75 K, and 
the stellar radius = 8 x 10 8 cm, so the model is determined only by the inner and outer 
radii of the disk. If we furthermore set the inner radius to be that at which dust sublimates, 
the only free parameter is the outer radius of the disk. Since emission from the inner radius 
dominates at the shortest wavelengths, the optically thick model makes robust predictions of 
the disk flux at the wavelength where emission from material at the vaporization temperature 
peaks, i.e. around 3 /im. The flux at longer wavelengths depends on the outer radius. 

As observed from Earth , the disk may of course be inclined with respect to the line of 



sight (Ivon Hippe 



illustration (IJura 



et al. 



2003, 



2007), though previous calculations considered face-on geometry for 



20071 ). The flux will scale as cosi, where i is the angle between the 
disk axis and the line of sight, until i —> tan^ 1 Hj Router, where Router is the outer radius of 
the disk an d H is its scale height. Taking the outer radius ~ 50i?* from the GD 362 model 
(I Jural 120071 ). and requiring disk thickness H < R*, the low-inclination limit only applies 



when i < 1.2°. The nearly edge-on limit is relatively improbable and furthermore would 
cover the star unless H << i?*, in which case the edge-on limit applies to even less probable 
geometries. An optically thick disk has a spectral energy distribution determined almost 
entirely by the outer radius, and total flux scaling with cosi. 

Using the Spitzer spectra, it is clear the simple optically thick model is definitively ruled 
out by the presence of a very strong silicate emission feature: this feature requires an optically 
thin emitting region. The feature contributes a significant portion of the disk luminosity and 
must have associated continuum. Further, the observed spectral energy distr ibution requires 
a colder component with color temperature ~ 290 K (jReach et all l2005al ) to explain the 
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flux at 24 /nil. As a first step toward constraining a possible optically thick disk around 
the star, we fitted the thin-disk model to the spectrum excluding the silicate feature (8-12 
jum), and setting the inner radius as the location where the grain temperature is 1200 K. 
The best fit has xl — 3-8; this high value is due to the lack of significant emergent cold flux 
from the model. However we take the constraints on R2 and % as a guide, with best values 
Ra/R. = 49 ± 5 and % = 41 ± 3°. 

Instead of attempting to fit the entire spectrum with the thin-disk model, we now 
consider only fitting it to the shorter-wavelength continuum, with an eye toward adding a 
cooler, physically-thick silicate-bearing cloud. Thus we excluded wavelengths longer than 
8 jum and fitted the physically thin model to the spectrum of G29-38. The near-infrared 
spectrum, and the decrease from 5-8 jum, require the inner edge of the disk is closer to the 
star than the point that reaches 1200 K; a better fit is obtained with T vap = 1500 K so 
that R\ = 9i?*. Figure [T2l shows the constraints on outer radius (R2) and inclination (z). 
The best fit has a low xl = 0-8; it is probably less than 1 due to overestimation of the 
uncertainties due to photospheric subtraction in the near-infrared. The best fitting outer 
radius R2 = 25-R* and inclination 23° . It is notable that the required inclination is roughly in 



the range required by iGraham et all (119901 ) in his model for the near- infrared timing (which 
required an inclination such that the dust temperature pulsations are detectable while the 
exciting pulsations on the photosphere are not). 

An amendment to the physically thin disk model is needed to improve the fit and 
explain the silicate emission. One possible solution is t o include an optically thin, effectively- 



physically-thick region at the outer edges of the ring. lJural (120071 ) showed that for GD 362, 
an extension of the ring that is warped, by only ~ 7°, can produce a silicate feature in 
emission. Figure [6] shows their model, compared to the moderately thick model from §4.1.21 
Indeed this outer region produces both the silicate feature and a significant fraction of the 
continuum at A > 11 fim in their model. The p lausibility of the warp of the outer disk is 



discussed and justified physically by I Jural (120071 ). who label it as 'region III' in their model. 



In many ways, the details (mass, temperature, radius, vertical extent) of this outer portion 
of the disk must be similar to the physically thick model, since they explain the mid-infrared 
emission in the same way (an optically thin cloud of silicates). In the warped disk model, 
the optically thin region is actually the upper layer of the disk; i.e. even the warped portion 
of the disk may be optically thick, as long as it is has a direct line of sight to the star. 
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4.3. Comparing the thick and thin models 

The disk spectral energy distribution can be empirically decomposed into three major 
components. One of them is a continuum with a relatively hot (890 K) color temperature, 
peaking around 4 /im and dominating the near-infrared emission. Another component is 
continuum with a lower color te mpera t ure ( ~ 300 K). The other major component is the 



silicate emission feature. In the lJural (120071 ) model, the hot component is the blackbody 



disk — hot because it is close to the star — and the cooler continuum and silicate feature arise in 
the outer warp region III. In the physically thick model (in Paper I and §4. 1.21 above), the hot 
component is amorphous carbon — hot because it is due to highly absorbing material — and 
the cooler component is silicates — cool because the silicates are more transparent and have 
strong mid-infrared emission features that allow them to cool efficiently; the two materials 
are colocated. 

For many plausible configurations of solid material around the white dwarf, we can 
consider the cloud as the sum of an optical ly thick disk and an optically thin halo or flared 



disk surface. Indeed, IVinkovic et all (120031 ) proved mathematically that flared disk models 
are equivalent to disk plus halo models. The radial profile of a spherical halo can be directly 
related to the flaring angle of a flared disk. Thus it is not possible, using the spectral energy 
distribution alone, to separate disk and halo (or warped disk) models. The halo dominates 
the infrared e mission when Thain > H/4R, where H is the flare height and R the distance 



from the star ( Vinkovic et al. 20031 ). For the wedge-shaped 'physically thick' model discussed 



above, this is equivalent to tii > tan#i/4. Equating the emission and absorbed flux yoyo 
of the stellar flux), the halo will dominate the infrared emission when tii > 0.004, which is 
already required for both the physically thick and thin models. 

For the specific case of the G29-38 disk, we made a direct substitution of the optically- 
thick disk for the C/Fe component of the moderately-optically-thick model. The optically- 
thick disk was taken directly from the model in Figure [TBI (i.e. the one fitted to wavelengths 
shorter than 8 /im), and the silicate components were taken from the moderately optically 
thick model in Figure [TDJ The silicate component of the moderately-optically-thick disk was 
rescaled in amplitude in order to match the observed flux after being added to the optically 
thick disk. The parameters for the best-fitting disk in this combined disk+silicate model are 
similar to those derived above using only short- wavelength data; R2 = 11R* and i = 29°. 
Figure dH shows the best fit; it is very similar to the moderately optically thick model in 
Figure [10| the goodness-of-fit, xl — 1-40, not much worse than the moderately optically 
thick model. Improvements to xl as obtained in §3] could be obtained by including more 
minerals in the optically thin region. 
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5. Timing constraints 



The geometry of the G29-38 disk can been constrained using timing info rmation. The 
star is a ZZ Ceti variable with non-radial modes that yield its optical pulsations. iGraham et al 
(Il990l ) found pulsations at periods ~ 200 sec in the K and L bands that have no strong coun- 
terpart in J and B-bands. More detailed optical photometry clearly shows pulsations at the 
predicte d frequency: see the peak around 5380 yuHz from the Whole Earth Telescope obser- 
vations (IWinget et al.lll990l ). The weakness of this mode at visible wavelengths, compared 
to infrared wavelengths, could be due to the manner in which it interacts with the disk. The 
B-band light is completely dominated by direct photons from the star, while the K-band light 
contains a contribution from thermal emission by the dust that produces the mid-infrared 
excess. The K-band pulsations are attributed to dust temperature variations. 



As I Graham et al.l (119901 ) explained, to generate observable K and L-band pulsations 
from a mode that is not prominent on the photosphere, the geometry is constrained. They 
invoke a mode of stellar brightness variations (2 nd order spherical harmonic) that brightens 
the pole and dims the equator of the star. If the disk is thin enough that its heating is 
driven by the equatorial stellar brightness, then it wi ll track the equatoria l photospheric 



temperature, rather than the average over the surface. I Graham et al. 



(119901 1 show that the 



spherical harmonic that excites the disk temperature variations can exist without significant 
variation of optical wavelength light if the star is viewed at the angle where the ratio of the 
two spherical harmonics is small. To eliminate J-band variability in their data required a 
viewing angle in the range 45-65°. 

The IRAC observations confirm the significance of the 200 sec period in the near- 
infrared. (No modes with periods longer than 500 s are detectable due to the brevity of 
the IRAC observ ations, so we could not address the 615 sec period that dominates the op- 
tical pulsations (IWinget et al.l Il990l ).) The fluctuating portion of the infrared flux, after 
subtracting the direct stellar photosphere, is 7% at 2.2 //m, 5% at 3.6 /im, and < 3% at 8 
fim. The color of the fluctuations suggests they arise from the inner portion of the disk, near 
the dust vaporization temperature. This explanation works well with the physically thin disk 
model with the inner edge of the disk at 1200 K and decreasing temperature outwards. The 
timing constraints rule out a spherical distribution for all the dust, though some dust can 
remain in a spherical distribution without generating pulsations (because the temperature 
variations average out over a sphere). The observations do not yet rule out a 'physically 
thick' configuration with a small opening angle (0i). We consider the IGraham et al.l (Il990l ) 
results very important; they require observational confirmation and have the promise of 
revealing the disk structure in more detail. 
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6. Discussion and Conclusions 

The mid- infrared spectrum of G29-38 is due to a cloud of small particles orbiting the 
star within the distance where tidal forces from white dwarf's strong gravity would destroy 
a large, weak body The tidal forces exceed self-gravity and strength of a non-rotating, rigid 
spherical body within a distance 



d/R* = 120p 



-1/3 



1 + 0.11^| 



(7) 



where p is the bodies density (g cm -3 ), r 10 is its radius (in units of 10 km), and S'kPa is its 
strength in kPa. If the dust around G29-38 was produced by one body, then the present- 
day mass requires r w > 1. For an asteroid or comet, we expect 0.5 < p < 2. Spin rates 

— 1/2 

of near-E arth objects req uire strengths ~ 2r 10 kPa to balance centrifugal forces against 



strength (jHolsappld 120071 ). The strength term in equation [7] is thus negligible for bodies 
with r w p > 1, which is likely for the parent bodies of the G29-38 dust. Fragments of parent 
bodies can survive much closer to the star; rocks with Sk pa = 3000, p = 3, and could survive 
against tidal disruption all the way to the surface of the star if they are smaller than 0.1 km. 

In the spectrum of G29-38, the strong emission feature at 9-11 /im, the color temper- 
ature, and the timing information all support the interpretation of the mid-infrared excess 
as being dominated by a cloud of small silicate particles. The infrared excess at 2-6 /im is 
apparently due to a featureless blackbody continuum. Based on our modeling results, this 
higher-color-temperature emission could be explained either by a highly-absorbing mineral 
(like solid C or Fe) or by a massive disk of material that is thinner than the white dwarf's 
diameter. The difference between these models, both of which can explain the observed 
spectral energy distribution rea sonably we ll, is significant: the massive disk model could 



harbor 10 24 g (or more) of dust (I Jural 120031 ) . while the cloud of amorphous C or Fe requires 
~ 10 19 g. There is no a priori reason that either of these mass estimates should be preferred 
or rejected. Some 10 19 g of silicates are required in all models. 

To understand the origin of the circumstellar material, it is important to know how 
much mass is observed. If only 10 19 g of material is present, then the observed infrared 
excess can be explained by a single small asteroid or a comet. If 10 24 g of material is present, 
then an entire, large asteroid, or numerous smaller ones, is required. In all cases, the parent 
body must have been somehow transported from a distance far enough from the star that it 
could have survived the red giant and AGB phase of the star (> 5 x 10 13 cm) inward to the 
Roche limit where it would be disrupted by the gravity of the white dwarf (< 10 11 cm). 
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6.1. Nature of the parent bodies 

The mineralogical results can help us relate the dust to possible parent bodies. We 
concentrate on the two primary components of the spectrum separately: the silicate material 
(required in all models in order to explain the 9-11 /zm emission feature) and the highly- 
absorbing material (required only in the thick-disk models). 

The composition of the highly-absorbing material cannot be determined unambiguously 
from the spectrum. We fitted it in $3]with amorphous C, based on the high cosmic abundance 
of C. But solid Fe or Si are also plausible, given their high cosmic abundance. Indeed, for 
chondrit ic (asteroid al and terrestrial) material the abundance of Si and Fe is much higher 
than C (jjura 2006 ). Mineralogical models for dust around other stars and in co mets, using 



the sa me methodology as in Sj3], do not always show a high abundance of C (ILisse et al. 



2007al lbh. For extrasolar systems, a careful subtraction of the photosphere is critical to 



measuring the high-color-temperature emission (which is characteristic of highly-absorbing 
material like solid C). For G29-38, the high-color-temperature excess (at 3.6 jum and longer 
wavelengths) is so far above photosphere that it must arise from circumstellar material, but 
at wavelengths shorter than 2 /im the infrared excess shape depends on the photosphere 
model. For G29-38 the mid-infrared data clearly require emission with a color temperature 
> 800 K, whether it is C or Fe. 

The compositions of the potential parent bodies for the circumstellar material around 
white dwarfs can be addressed by studies of Solar System bodies. To date there has been 
no sample return mission from an asteroid, but meteorites provide direct measures of com- 
position of parts of at least some asteroids. Carbonaceous chondrites have some C, but 
all chondrites are largely silicate mineral (olivine and pyroxene, mostly Mg-rich [forsterite 
and enstatite]), with a wide range of other minerals (some Ca and Al-rich) and metals (of- 
ten including previously- molten Fe). Metallic meteo rites, commonly found in museums an d 



on the ground, have largely Fe and Ni composition (IShearer. Papike. and Rietmeijerlll998l ). 
Comets are likely to have a more primitive composition than carbonaceous chondrites, with 
abundant silicate grains as well as carbonaceous material, based on infrared spectroscopy of 
cometary dust, laboratory study of cometary inte rplanetary dust particles , and in situ mass 



spectrometry during the lP/Halley flyby in 1986 (lHanner fe Bradleyll2004l ). The most abun- 
dant silicate minerals in meteorites are Mg-rich olivines and pyroxenes, as well as feldspar 
and phyllosilicates. Interplanetary dust particles believed to originate from comets (CP type) 
are largely composed of phyllosilicates th at require aqueous alteration on their parent body 



([Messenger. Sandford. fe Brownleel 120061 ) . 



Based on the analogy to Solar System bodies, the dominance of Fe-rich pyroxene min- 
eralogy is distinct. Comets or outer-main-belt (D-type) asteroids contain Mg-rich olivine 
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and phyllosilicates, which when combined can reasonably fit the observed spectrum of G29- 
38. Both comets and D-type asteroids contain organic material, which would be consistent 
with the presence of amorphous C, and the most primitive carbonaceous (CI) meteorites are 
largely composed of phyllosilicates. But the Mg-rich olivine plus phyllosilicate model is not 
as good a fit to the data as Fe-rich pyroxene. There are pyroxene-dominated meteorites, but 
they are dominated by Mg-rich pyroxene (hence the name enstatite meteorites for this rare 
class). It is worth noting that enstatite chondrites are thoug ht to have formed in r educing 
conditions and contain other minerals including niningerite (jWeisberg et al.ll2006l ). which 
was one of the most abundant minerals from our fit for G29-38 in Tab. [31 Thus the physical 
conditions for formation of the enstatite meteorite parent bodies may have some relevance 
to the formation of the G29-38 debris parent body. 

In terms of the featureless material that produces the near-infrared continuum (in the 
physically thick models), either C or Fe are acceptable to the fits. If the material were C, 
then a cometary or D-type asteroid origin would be more likely, whereas with abundant Fe, 
formation closer to the star and within a differentiated parent body would be implicated. 

Some constraints on the composition of the material are obtained from the abundances 
in the white dwarf atmosphere. Solid material at the inner edge of the disk is constantly 
being vaporized by the stellar radiation. These vapors reach the atmosphere of the star 
(or are blown out of the system). They cannot reside in the stellar atmosphere for long; 
instead, they diffuse rapidly inward, deeper than the photosphere. Thus the heavy elements 
in the stellar atmosphere must be 'fresh,' consistent with an origin from vaporization of 
circumstellar dust but not wi th a long-lived stellar atmosphere. G29-38 has metals present 
in its atmosphere. lJural (120061 ) discussed the deficiency of C in some externally polluted white 
dwarfs, where the abundance of C relative to Fe is more than 10 times lower than solar. CI 
meteorites have C/Fe 10 times lower than solar; this is commonly explained by th e volatility 
of C and the high temperature required for chondrite formation (lBrearleylll998l ). We note 
that the abundances of refractory elements in comets and asteroids as w ell as that inferre d 
from exozodiacal dust is consistently less than solar, in the 7-10 % range (ILisse et al.ll2006l ). 
But in a relative abundance, dust collected during the encounter with the lon g-period comet 
lP/Halley has C/Fe abundance ratio similar to Solar (IJessberger et al.lll988l ). Short-period 
comets (which are periodically heated to higher temperature) may be more devolatilized and 
may have abundances more similar to asteroidal material. 
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6.2. Disruption of 'Hot Jupiter'? 



One intriguing possibility for the origin of the infrared excess is the survival of the 
core of a giant planet and its subsequent gravitational disruption. 'Hot Jupiters,' with 
masses of order 10 30 g orbiting within 10 12 cm (0.1 AU) of their star, appear to be fairly 
common in extraso lar planetary systems: 1.2% of nearby F, G, and K stars has such a planet 
jMarcv et al.ll2005h . In this scenario, the planet would become engulfed into a common 
atmosphere during the red giant phase. Drag from the extended stellar atmosphere would 
cause the planet to spiral inward toward the stellar core. At the end of the mass-losing phase 
of the star's evolution, we would be left with a white dwarf, the surviving core of the planet, 
and the planetary nebula composed of the outer atmosphere from the star. 

If a planet began at < 10 12 cm from the star, it is possible that its remnant could land 
within the Roche radius. The effect of the post-main sequence evolution on the planetary 
dynamics has not been explored in detail. The change in mas s of the central star could m ake 
any borderline- unstable system of multiple planets unstable ( IDebes fc Sigurdssonl 120021 ) and 
could lead to nonlinear orbital perturbations. A simple estimate of the orbital decay due to 
gas drag when the planet is within the red giant atmosphere is made by setting the rate of 
kinetic energy imparted to the planet, 

3/2 



E = l - P TtR 2 v 
equal to the change in orbital binding energy 



GM* 



U 



GM*M P da 
a~ 2 d? 



(9) 



where M p and R p are the mass and radius of the planet, G is the gravitational constant, a is 
semimajor axis of the presumed-circular orbit, M* is the mass of the star, and p is the mass 
density of the star at the distance of the planet. Taking for illustration a planet with the 
mass of Jupiter and average density 1 g cm -3 , orbiting at a distance of 0.1 AU from a star 
of mass 2M that has mass density 10 -9 g cm -3 at 0.1 AU from its center, the timescale for 
orbital dec ay is 6 x 10 4 yr. This is mu ch smaller than the asymptotic giant br anch lifetime 
10 7 yr (jVassiliadis and Wood! Il993l ) for a progenitor star mass of 3.1M Q (IWeidemann 



20001 ) . Thus the inward spiral of hot Jupiter (or other close-in) planetary orbits is likely to 
be rapid, once the star has entered its giant phase. 



Considering only a convective main-sequence primary star, ISasselovl (120031 ) showed that 
the orbit of a 'hot Jupiter' like OGLE-TR-56b would tidally decay on a timescale of 1-10 
Gyr. Again considering only a main sequence primary star, iBaraffe et al.l (120041 ) showed 
that planets with masses less than 3 x 10 30 will evaporate down to a rocky core in less than 
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5 Gy r These effects may combine, with tidal decay enhancing evaporation (jErkaev et al. 
20071 ). Tidal decay of 'hot Jupiter' orbits cannot be too rapid for main sequence stars, or 
else they would be much less common. The orbital evolution is very sensitive to the stellar 
structure and the planet's orbital eccentricity; the response is highly nonlinear and may lead 
to a rapid inward spiral upon perturbation of the stellar interior structure or eccentricity 
( lOgilvie and Linl 120071 ; IJackson et al.l 120071 ) Nonetheless, regardless of whether tidal decay 
can bring planets close to the stellar surface during main sequence evolution, the gas drag 
during the red giant phase should finish the job. 

How much of a 'hot jupiter' could survive the late stages of its star's evolution is not 
known. A brown dwarf was recently discovered around a white dwarf suggesting that some 
companions can survive the red g iant phase of the primary star despite being well within the 
atmosphere (IMaxted et al.ll2006l ). Indeed, the companion is within the Roche distance of the 
white dwarf unless it has a mass density greater than 3.6 g cm -3 or is held together by more 
than its own gravity. The composition of the core of an extrasolar giant planet is difficult to 
constrain. Mass and radiu s estimates seem to r equire a dense core, possibly of high-pressure 
ice or stony composition (jBurrows et al.. 1 120071 ). This material would be tidally disrupted, 
as in the asteroid (or comet) disruption models considered to date. The Fe-rich pyroxene 
mineralogy for the G29-38 debris, which is distinct from that of most comets and asteroids, 
does not appear impossible for a rocky planetary core, but at present there is little more 
than can be said about whether such a mineralogy is likely. 

That G29-38 is exceptional among debris disks (having by far the brightest infrared 
excess, with 3% of the star's luminosity absorbed and radiated by dust) makes the hypoth- 
esis more attractive. Consider the reverse argument: if 1% of all stars have 'hot jupiters,' 
what is their expected end state? From this point of view, it seems inevitable that such 
remnant planets would generate debris for their parent star in its white dwarf phase. But 
this argument assumes that 'hot jupiters' exist around A-type stars (the progenitor type for 
G29-38), at least part of the planet can survive the common envelope phase, the remnant 
core lands within the Roche radius, and the debris can survive 500 Myr (the age of G29-38). 
The age problem may not be severe, if the remnant planetary core becomes fragmented, with 
fragments gradually entering the Roche radius due to collisional disruption and gravitational 
perturbations. Further theoretical work is needed to test the viability of this scenario. 
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Table 1: Spitzer Observing log for G29-38 

Date AORID Instrument Wavelengths 

2004 Nov 26 10119424 IRAC 4.5, 8 /xm 

2005 Nov 26 11124224 IRAC 3.6, 4.5, 5.8, 8 /im 
2004 Dec 2 10149376 MIPS 24 fim 

2004 Dec 8 10184192 IRS 5.2-14.2 /im 

2005 Dec 23 13835264 IRAC 3.6, 4.5, 5.8, 8 fim 

2006 Jun 30 13828096 IRS 5.2-36 /iin 
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Table 2: Mid-infrared Fluxes of G29-39 (mJy) a 

2004 Nov 26 2004 Nov 26 2005 Dec 23 
10:54 10:58 23:23 

3.6 /im 8.37 ±0.01 8.10 ± 0.03 

4.5 /im 8.88 ±0.02 8.87 ±0.01 
5.8 //m 8.37 ±0.02 8.28 ± 0.07 

8.0 //m 8.73 ±0.03 8.72 ± 0.02 
"Uncertainties are statistical uncertainty in the weighted mean of the flux measurements 
from all frames taken during each observing sequence 
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Table 3. Composition of the Best-Fit Model a to the Spitzer/IRS G29-38 Spectrum 



Species 


Weighted 15 


Density 


Mol. Wt. 




± max 


X 2 if 




Surface Area 


(g cm" 3 ) 


(relative) 




(K) 


excluded 






Del 


sections 










Amorph Olivine (MgFeSi0 4 ) 


0.33 


3.6 


172 


0.69 


890 


90.6 


Fayalite (Fe 2 Si0 4 ) 


0.08 


4.3 


204 


0.17 


890 


2.91 


FerroSilite (Fe 2 Si20 6 ) 


0.11 


4.0 


264 


0.17 


890 


9.85 


Diopside (CaMgSi 2 6 ) 


0.05 


3.3 


216 


0.076 


890 


2.05 


OrthoEnstatite (Mg 2 Si 2 6 ) 


0.04 


3.2 


200 


0.064 


890 


1.98 


Niningerite (Mgi FegoS) c 


0.10 


4.5 


84 


0.53 


890 


1.49 


Amorph Carbon (C) 


0.28 


2.5 


12 


5.83 


930 


> 100 


Water-ice (H 2 0) 


0.29 


1.0 


18 


1.61 


220 


5.82 


Upper Limits and Non-Detections 


Forsterite [Koike] (Mg 2 Si0 4 ) 


0.02 


3.2 


140 


0.046 


890 


1.15 


Amorph Pyroxene (MgFeSi 2 Og) 


0.00 


3.5 


232 


0.09 


890 


1.04 


Smectite/Notronite' 


0.00 


2.3 


496 


0.03 


890 


1.04 


Water Gas (H 2 0) 


0.01 


1.0 


18 


< 0.00 


220 


1.04 


Magnesite (MgC0 3 ) 


0.00 


3.1 


84 


< 0.00 


890 


1.04 


Siderite (FeC0 3 ) 


0.00 


3.9 


116 


< 0.00 


890 


1.04 


PAH (C10H14) 


0.00 


1.0 


(178) 


< 0.011 


N/A 


1.04 



a Best-fit model xl = 1-04 with power law particle size distribution dn/da oc a 3 7 , 5-35 /im range of fit, 
336 degrees of freedom 

b Weight of the emissivity spectrum of each dust species required to match the G29-38 emissivity spectrum. 

c N mo i es (i) is the Density/Molecular Weight x Normalized Surface Area for mineral i. Errors arc ±15% 
(la). 

d All temperatures are ±20K (la) 

e We use the name niningerite to refer to Mg^Fei-^S. a niningerite composition of Mg 25 Fe75S may fit the 
data better. 

f Na .33Fe 2 (Si,Al)4Oi (OH) 2 . 3 h 2 0) 
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Table 4: Best-fitting moderately-optically-thin model 

Parameter Best value Confidence interval 

1050 - 1200 K 
1 - 8 
< 150 
2.4 - 2.9 



T vap 1100 K 

T|| 2 

R-min 50 

a 2.7 
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Fig. 1. — Flux versus time for each of the IRAC channels (from top to bottom: 3.6, 4.5, 5.8, 
and 8 /zm wavelength). The time begins with the first frame of the entire sequence (which 
combines two AORs) at 2004 Nov 26 10:54:11.28 UT. Each point is a flux measurement 
from a single image, with the statistical uncertai nty of the aperture photometry performed 



as described by the IRAC calibration procedure (IReach et al. 



2005bl ). Line segments simply 



connect the data points. The high-error points at 4.5 and 8 /im are the short frames taken 
just at the beginning of each of the two AORs. All four plots are scaled to show the same 
dynamic range (from -6% to +8.4%, centered on the median). 
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Fig. 2. — Scargle periodograms of the time series from the IRAC 2004 Nov 26 observations. 
Each panel has the same scale both horizontally and vertically; the 3.6 fxm panel is twice as 
large because the amplitudes are much higher. The periods (2ir/u, where u is the angular 
frequency) are displayed on a logarithmic stretch. 



-36- 




5 10 15 20 25 50 
wavelength (/xm) 



Fig. 3.— Combined spectra and photometry (using 2 MASS, IRTF, IRAC, MIPS, and IRS) 
of G29-38 with stellar photosphere removed. 
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Fig. 4. — Emissivity spectrum of the dust around G29-38. The observed spectrum has been 
divided by a 930 K blackbody and fitted with a linear combination of 12 minerals. The 
shape of each mineral's emissivity, normalized by its fitted amplitude to the G29-38 emis- 
sivity, is shown separately by a colored line (offset vertically for clarity): red = amorphous 
carbon, bright green = carbonates (zero amplitude), yellow = PAH (zero amplitude), light 
orange = water vapor (zero amplitude), deep orange = water ice, olive green = sulfides, rep- 
resented here by niningerite, blue = phyllosilicates (zero amplitude), light blue = crystalline 
pyroxenes (ferrosilite, diopside, and orthoenstatite, in order of 20 /im amplitude), purple = 
amorphous olivine, and dark blue = crystalline olivines (forsterite and fayalite, in order of 
20 /im amplitude). 



-38- 




5 10 15 20 25 30 35 

Wavelength (^im) 

Fig. 5. — Emissivity spectrum as in Fig. HI but after removal of the best-fitting silicates and 
carbon. The residuals are well fit by a combination of water ice (deep orange, dashed and 
solid are two temperatures) and metal sulfides (olive green). The water ice is at temperature 
200 K, and cannot be in the same location as the 930 K dust. Whereas the dust is ~ 10 11 
cm from the star, the water ice must be further ~ 10 13 cm, at the outer edge of an extended 
disk. 
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Fig. 6. — Cartoons illustrating possi ble geometr ies for the G29-38 disk. The physically 
thin plus warp model (top) is based on lJural (120071) 's model for the s imilar white dwarf GD 
362. The physically thick (bottom) is based on iReach et al.l (j2005al ) and elaborated in §4j 
Distances from the star are labeled, in units of stellar radii. The distances from the star for 
grain vaporization (greybody grains at 1200 K) and tidal disruption (Roche limit for a solid 
body to be tidally disrupted by the white dwarf with radius 7.5 x 10 8 cm and mass 0.69 M Q ) 
are indicated. 
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Fig. 7. — Spherical shell models for the dust around G29-38. Individual models for amor- 
phous olivine and carbon are shown as dashed and dotted curves, respectively, and a lin- 
ear combination is shown as the solid curve. The input white dwarf spectrum, and the 
photosphere-subtracted infrared observations, are shown for comparison. The DUSTY mod- 
els have optical depth r(0.55/xm) = 0.018 and 0.011 in silicates and carbon, respectively. 



-41 - 




Fig. 8. — Temperature profiles in a spherical shell around G29-38. The top panel is for 
amorphous olivine, and the bottom panel is for amorphous carbon. All models are computed 
to an inner temperature of 1200 K. Each profile is for a different total optical depth, with 
values 0.0100, 0.0268, 0.072, 0.193, 0.518, 1.39, 3.73, and 10. The optically thin models 
extend closest to the star, are the warmest at the outer edge of the plot, and are closest 
to a straight line (T oc r~ 0,5 ). The optically thick models begin further from the star and 
decrease in temperature much more steeply. 
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Fig. 9. — Temperatures of olivine (top) and carbon (bottom) grains of three different sizes: 
0.1 /im (diamonds), 1 /im (triangles), and 10 /zm (squares). The grains are irradiated directly 
by the white dwarf (i.e. the cloud is optically thin). The solid line shows power-law fits of the 
form T = Ti(r/100i?*)~ 5 . For olivine (carbon) grains of 0.1, 1, and 10 /im radius, 7\ = 1905, 
1136, 853 K (1393, 740, 780 K) and 5 = 0.47, 0.50, 0.54 (0.37, 0.41, 0.51), respectively. The 
horizontal dotted line indicates the vaporization temperature (1200 K) in our model. 



-43- 




10 20 30 40 

Wavelength (/um) 



Fig. 10. — Best-fitting moderately-optically thick disk model (solid line) combining amor- 
phous carbon (dashed), amorphous olivine (dotted), crystalline bronzite (dash-dot) to the 
observed photosphere-subtracted spectral energy distribution of G29-38. 
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Fig. 11. — Residuals (observed minus model) from the best-fitting moderately-optically 
thick disk model of Fig. [TUJ Features discussed in the text are labeled. 




Fig. 12. — Goodness of fit (x 2 /d.o.f) from the physically-thin, optically-thick model for a 
range of outer radii (R2, in units of stellar radius) and inclination. The best fitting model 
has inner and outer radii R\ — 9, R 2 — 25, and inclination 23°. 
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Fig. 13. — Best-fitting physically-thin, optically-thick disk model. The model (solid line) 
was only fitted to data at wavelengths less than 8 /xm; it cannot produce the silicate feature 
and under-predicts (by far) the longer-wavelength continuum. 
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Fig. 14. — Best-fitting combined model, with a physically-thin, optically-thick disk (cf. 
Fig [13]) plus a silicate-only (no C or Fe) physically-thick disk (cf. Fig. [TOl) . 



